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Effective Dynamics for the Cosmological Bounces 
in Bianchi Type I Loop Quantum Cosmology 

Dah-Wei Chiou 

Institute for Gravitational Physics and Geometry, Physics Department, 
The Pennsylvania State University, University Park, PA 16802, U.S.A. 

The detailed formulation for loop quantum cosmology (LQC) in the Bianchi I model with a 
massless scalar field was recently constructed. In this paper, its effective dynamics with the LQC 
discreteness corrections is studied and the equations of motion are analytically solved, showing that 
the big bang singularity is replaced by the big bounces, which take place up to three times, once in 
each diagonal direction, whenever each of the area scale factors approaches its critical value in the 
Planck regime measured by the reference of the scalar field momentum. 
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I. INTRODUCTION 

The comprehensive formulation for loop quantum cos- 
mology (LQC) in the spatially flat-isotropic model has 
been constructed [H, [|J • With a massless scalar field serv- 
ing as the emergent time, the result shows that the quan- 
tum evolution is deterministic across the deep Planck 
regime and in the backward evolution of the states which 
are semiclassical at late times, the big bang is replaced by 
a big bounce. Based on the same principles, the construc- 
tion was further improved by a more direct implementa- 
tion of the underlying physical ideas of loop quantum 
gravity (LQG) [3]. In the improved dynamics, the big 
bounce occurs precisely when the matter density enters 
the Planck regime, regardless of the value of the momen- 
tum p,), of the scalar field. 

Both the precursor strategy ( "/i -scheme" ) and the im- 
proved strategy ("^-scheme") were applied and recon- 
structed for the Bianchi I model to include anisotropy 
The analytical investigation shows that the state in 
the kinematical Hilbert space associated with the clas- 
sical singularity is completely decoupled in the difference 
evolution equation, indicating that the classical singu- 
larity is resolved in the quantum evolution and the big 
bounce may take place when any of the area scales un- 
dergoes the vanishing behavior. 

While a thorough numerical investigation remains to 
be done to draw the definite conclusion for the details of 
the quantum evolution in the Bianchi I model, this pa- 
per studies its effective dynamics with LQC discreteness 
corrections. Not only does the result affirm the anticipa- 
tions in 0] but more intuitive pictures are also obtained 
in this semiclassical approach, giving an insight into how 
and why the big bounces take place. 

In accordance with the formulation in Q , this paper fo- 
cus specifically on the model with a massless scalar field. 
In the context of effective dynamics with LQC discrete- 
ness corrections, the similar analysis for more generic 
Bianchi I models with the inclusion of arbitrary mat- 
ter with equation of state w < +1 is also investigated 
in H, which gives the similar results for the occurrence 
of big bounces with only difference in detail. With arbi- 



trary matter sources, however, the equations of motion 
are very complicated and a proper approximation has to 
be used. By contrast, in the special case of a massless 
scalar field, the equations of motion can be solved an- 
alytically and therefore the underlying physics is more 
transparent. 

This paper is organized as follows. In Section HH the 
classical dynamics of the Bianchi I cosmology with a 
massless scalar source is solved in terms of Ashtekar vari- 
ables in Hamiltonian formulation. The effective dynam- 
ics with LQC corrections in /2-scheme is constructed and 
solved in Section IIIII Its phenomenological ramifications 
are discussed in Section IIV1 As a comparison to the p- 
scheme, the effective dynamics in // Q -scheme is also in- 
cluded in Appendix [Al 



II. CLASSICAL DYNAMICS 

The spacetime metric of Bianchi type I is given as: 

ds 2 = -dt 2 + a\(t)dx 2 + a 2 2 (t)dy 2 + al(t)dz 2 . (2.1) 

In terms of Ashtekar variables, the phase space of Bianchi 
I models is given by the diagonal triad variables pi and 
diagonal connection variables cj for / = 1,2,3, which 
satisfy the canonical relation: 



{ci,pj} = 87rG7<5/j. 



(2.2) 



The triad variables pj are related with the length scale 
factors aj via: 



Pi = a 2 a 3 , 



P2 = aia 3 , 



P3 = aia-2- 



(2.3) 



In the presence of a massless scalar field <j>(x,t) = 4>(t), 
(which is independent of the spatial coordinates with ho- 
mogeneity assumed) , the classical dynamics is govern by 
the Hamiltonian constraint: 



C — Cg ra v + C(j) 



(2.4) 



(c 2 p 2 C 3 P3 + CiPiC 3 p 3 + C1P1C2P2) 
%TTG~f 2 ^Pip 2 p 3 



Pj, 



2 



where p^, is the conjugate momentum of <f> and has the 
canonical relation with 6: 



{0,P4>} = 1. 



(2.5) 



We can simplify the Hamiltonian by choosing the lapse 
function N = y / piP2P3 and thus introducing the new 
time variable dt' = (j>\p 2 p3)~ 1 l 2 dt. The rescaled Hamil- 
tonian constraint is given by 



H 



(C2P2C3P3 + CiPiC 3 p 3 + C 1 p 1 C 2 p2) 

87rG 7 2 



The equations of motion are governed by the Hamil- 
ton's equations: 



dp4> 
~dY 

dt 
dc 
~db 



{p<p,H}=0 



- = {<P,H} = p < j>, 
dci r m s n dH 

= -l~ 1 Cl(c 2 P2+C 3 p 3 ) , 

"<" 1 i //i a n dH 
— = { Pl ,H} = -8irG~f — 

1~ Y V\ {c 2 p 2 + c 3 p 3 ) , 



is constant (2.7) 
(2.8) 

(2.9) 



dp 



(2.10) 



and so on for c%, c 3 , p 2l p 3 in the cyclic manner. In 
addition to the Hamilton's equations, the constraint that 
the Hamiltonian must vanish yields 



H(ci, Pl ) = => 



(2.11) 



1 



■{c 2 p 2 c 3 p 3 + CiPiC 3 p 3 + c 1 p 1 c 2 p 2 ) 



e * 47rG 7 2 
Combining (pH)]) and (f5TTU]l gives 

— (pid) = 0, p I c I =MG 1 hK. I , (2.12) 

where ICi are dimensionless constants, which will be used 
to parameterize the solutions of evolution. Taking (|2.12j) 
into (|2.1ip . we have 

p\ = l&vGh 2 {/C 2 /C 3 + /C1/C3 + K1K2} (2.13) 

or equivalently 



= 2 (/C2/C3 + AC1AC3 + AC1AC2) , (2.14) 



if we define 



p,p := hv8nG JC 4 
Putting (|2"7T2"|) into (|2~TU|) gives 



(2.15) 



i^l=87rGft(/C 2 +/C 3 ) 
pi dt' 



(2.16) 

By referring to (|2.8[) , this leads to 

l^L = 87 :Gh'^±^ = V^G( 1 —^), (2.17) 



where we scale the parameters ICi = JCki, K,^ 
such that 



K,K, d 



Ki + k 2 + k 3 = 1, k\ + k\ + kI + kI= 1. (2.18) 

Regarding <fi as the emergent time, the solutions of evo- 
lution are given by 



Pi{<t>) =Pi(4>o)e- 



/8ttG{ 



)(<!>-M 



_l t_ (2.6) or equivalently 



ai{4>) = ai(4>o) e 



/8¥G ^-{4>-<t>o) 



(2.19) 



(2.20) 



The classical Bianchi I model with a massless scalar 
field admits both "Kasner-like" (two of K] positive and 
the other negative) and "Kasner-unlike" (all K] positive) 
solutions. The Kasner-like solution, which has two ex- 
panding and one contracting directions (say > 0), 
eventually encounters the "Kasner-like singularity" (a 
given regular cubical cell stretches as an infinitely long 
line) in the far past and the "planar collapse" (a regu- 
lar cubical cell stretches as an infinitely large plane) in 
the far future. On the other hand, the Kasner-unlike 
solution, with all directions expanding, encounters the 
"Kasner-unlike singularity" (a regular cubical cell van- 
ishes to a point) in the far past and no planar collapse. 

We will see that with LQC discreteness corrections, 
both Kasner-like and Kasner-unlike singularities are re- 
solved and replaced by the big bounces, whereas the pla- 
nar collapse remains its destiny even one of the three di- 
agonal directions approaches infinitely small length scale. 



III. EFFECTIVE DYNAMICS IN ^-SCHEME 

In LQC, the connection variables cj do not exist and 
should be replace by holonomies. In the effective the- 
ory, to capture the quantum corrections, following the 
procedures used in the isotropic case 0, 0] , we take the 
prescription to replace cj by sin(/i/c/)//2/, introducing 
discreteness variables /2/. In the improved strategy (/2- 
scheme) used in Bianchi I LQC [|| , fij are not fixed con- 
stants but given by 



A*3 



where A = -^(47T7£pj) is the area gap in the full theory 
of LQC. 

Imposing this prescription plus the loop quantum cor- 
rection to the inverse triad on (|2.4j) , we have the effective 
Hamiltonian constraint to the leading order: 





C' 



/(pi)/(pa)/(ps) 



Pi f(Pl)f(P2)f( P3 ) 



Pi 



2 8ttG7 2 
sm(p, 2 c 2 ) singes) 



(3.2) 



P2V3 



P2P3 + cyclic terms 



3 



where f(pi) is the eigenvalue of the inverse triad oper- 
ator 1/yfpi. The loop quantization gives the quantum 



and 



corrections: 



f(Pi) 



Hp,(ci,Pi) = 



(3.9) 



oc P/74? +1 for Pl < 



>2 
PI 
>2 
PI 



(3.3) 



47rG 7 2 



sin(^ 2 c 2 ) sin(/x 3 c 3 ) 



P2P3 + cyclic terms 



with the Planck length £p\ :— \/GH and a positive n. The 
corrections to f(pi) is significant only in the Planckian 
region in the vicinity of pj — 0. From now on, we will 
ignore the quantum corrections to f(pi) by simply taking 
its classical function f(pi) — Pj ■ [We will see that in 
the backward evolution the big bounce takes place much 
earlier before the discreteness correction on the inverse 
triad operator becomes considerable, and it is the "non- 
locality" effect (i.e., using the holonomies) that accounts 
for the occurrence of the big bounce.] 

With f{p I ) = pj 1/2 , by choosing dt' = {p 1 p 2 pzY xl2 dt 1 
the Hamiltonian constraint (13.21) can be rescaled as 



[Note that in the classical limit /2/Cj — > 0, we have 
sm(fizCi) / fti cj, cos(/2/Cj) — > 1 and therefore (|3.7p ~ 
reduce to their classical counterparts J2.9|) - (|2.11D .] 
By (|3~7)) and (|3~5| . we have 



3sin(/ijcr) c/cos(u/C/)\ dpi ,_ s dcj 

dF +PlC0S ^ lCl) W 



2/2/ 



d_ 

IF 



vi- 



sin(^/cj) 



/./,/ 



= 0, 



which gives 



sin(^/c/) 

= 8n&jnK,i. 



(3.10) 



(3.11) 



Hil = Y 



(3.4) 



1 fsin(//2C2)sin(/i3C 3 ) 



87rG 7 2 



^2^3 



P2P3 + cyclic terms 



Again, the equations of motion are given by the Hamil- 
ton's equations and the constraint that the Hamiltonian 
must vanish: 



~dF 

d(f> 

dt> 



Taking (|3.11|) into (|3.9j) again gives the same constraints 
on the constant parameters as in (|2.13p or (|2.14p . 
Substituting ([3~TT]) into (|3~5|) yields 

— = 8nGhcos(ji 1 c 1 )(>C 2 +JC 3 ). (3.12) 
Pi dt' 

By regarding (f> as the emergent time via p.6[) and ex- 



Hp,} =pj, 



is constant (3.5) pressing cosx = ±y/l- sin 2 x, (|3.12|) then gives 



(3.6) 



1 tfoj 



= ±V87rG 



1 — Kl 

K4, 



1 - 



Pi 



-,1/2 



Ql, crit 



"7 



_! /3 sin(/iici) ciCOs(/iiCi) 



2/2i 2 

/ sin(/2 2 C2) , sin(/2 3 c 3 ) 

= Pi + — 

V M2 M3 



P3 , (3.7) 



where we define the directional density: 

Qi ■= -f 

Pi 



(3.13) 



(3.14) 



dpi r it \ c n dH » 

FF = = -fc 

= 7 _ Vi cos(/2iCi) 

/ sin(/2 2 C2) , sin(/2 3 c 3 ) 



for the /-direction and its critical value is given by the 
Planckian matter density pp\ times a numerical factor: 



V M2 



-P2 



M3 



glcrit : — I 

Kl 



PPh 



PP1 



SttG^A)- 1 . (3.15) 



P8 , (3-8) 



IV. DISCUSSION 

As opposed to the classical equation (|2.17[) . in which pi 
continues to decrease toward the classical singularity in 
the backward evolution, the effective equation in (|3.13p 
flips sign exactly at the moment when gi approaches its 
critical value £?/, C rit. Note that by (|2.3|1 pi can be re- 
garded as the area scale factors. Therefore, with the 



LQC discreteness corrections, (I3.13P shows that the sin- 
gularities (both Kasner-like and Kasner-unlike) are re- 
solved and replaced by the big bounces in the backward 
evolution when any of the area scales undergoes the van- 
ishing behavior. Across the bounces, the equation of mo- 
tion again comes closer and closer to the classical coun- 
terpart. Hence, the semiclassicality is retained on both 
asymptotic sides of the evolution. 
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Furthermore, the detailed evolutions of pi are decou- 
pled in different diagonal directions and evolve indepen- 
dently of one another once the initial conditions (pi(<f>o), 
Pfi and Kj) are specified. Thus, the bounces occur up 
to three times, once in each direction, whenever each of 
the directional densities gi approaches its critical value. 
As expected, in /^-scheme, the critical values gj^ cr jt are 
in the Planck regime of 0(Mp i ) and independent of the 
value of p,], (p^crit depend on p$ only through the ratio 
k$/ki = Ktjj/Ki). 1 Note that gi have the same dimen- 
sion as the matter density p := p J (2pip2P3) and gi play 
the same role as p does in the isotropic case, signaling 
the occurrence of big bounces. 

On the other hand, the planar collapse is not re- 
solved but one of the length scale factors ai continues 
the vanishing behavior in the Kasner-like case. This is 



1 In Appendix fS] the old precursor strategy (^t -scheme) is pre- 
sented and it shows that the critical value of qi can be made 
arbitrarily small by increasing . 



expected since the classical solutions (|2.12j) and (I2.19P 
yield pici — > (and pfci — > in /^ -scheme) toward the 
planar collapse and therefore the quantum corrections 
become more and more negligible (in both schemes). 

For given initial conditions, the differential equation 
P-13P can be solved numerically. The behaviors of Pi{4>), 
ai(<f>) and gi(4>) are depicted in parts (a), (b) and (c) 
respectively in FIG. [1] and FIG. [2] for Kasner- unlike solu- 
tions and in FIG. [3] and FIG. 2] for Kasner-like solutions. 

The fact that smallness of pi (not of aj) is an indica- 
tion of the occurrence of big bounces seems to support the 
suggestion that "area is more fundamental than length 
in LQG" , although whether this is simply a technical ar- 
tifact or reflects some deep physics is still not clear. (See 
Section VII. B of [8[ for some comments on this aspect 
and [9j for more details.) Meanwhile, as the length op- 
erator has been shown to have a discrete spectrum [lfj, 
the fact that the vanishing of the length scale factor in 
the planar collapse is not stopped seems to contradict 
the discreteness of the length spectrum. Whether we 
miss some important ingredients when imposing the fun- 
damental discreteness of LQG in the LQC construction 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

#-*„ (G- 1/2 ) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

*-*o (G- 1/2 ) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

*-*„ (G" 1/2 ) 



FIG 



3: «i = 1/2, k 2 = 3/4, K 3 = -1/4, «^ = l/VI; pi(<M = P2(0o) = Ps(^o) = 1 X lO 4 ^; and p^ = 2x 10 3 ft\/8^G (i.e., 



2 x 1(T 




(b) 500 



0.2 0.4 0.6 0.8 1 



if 200 




*o (G ' ) 



0.2 0.4 0.6 0.8 1 

*-*o (G" 1/2 ) 




0.2 0.4 0.6 0.8 1 

*-*„ (G" 1/2 ) 



FIG. 4: ki = 1/2, k 2 = 3/4, k 3 = -1/4, ^ = l/x/8; pi(0 o ) = 3 x 10 4 4i> PaOa) = 2 x 10 4 4i, P3(0o) = 1 x 10 4 4i5 and 
p = 2 x 10 3 fiVc^G (i.e., K.H4, = 2 x 10 3 ). 



or indeed area is more essential than length remains an 
open question for further investigation. 

It is also noteworthy that (|3.13[) remains invariant if we 
rescale p^ — > Z 3 p0 and pi — > Z 2 pj at the same time. This 
is reminiscent of the idea as suggested in that area 

is measurable only if the surface is coupled with the ma- 
terial reference. The scaling invariance, however, breaks 
down in the full LQC theory since the quantum evolu- 
tion is governed by a difference equation [J| , in which the 
step size of difference introduces an additional scale in 
the deep Planck regime. 2 

Meanwhile, related to the above observation, the phys- 
ical meaning of the directional densities qi can be inter- 
preted as the (inverse of) area scales, again, measured 
by the reference of the matter content. The big bounces 
take place whenever one of the area scales becomes very 



2 Therefore, the scmiclassicality is retained in the full quantum 
theory only for large pj, and pi- Accordingly, we put big values 
of p^ and pi(<fio) in the figures to make sense of the semiclassical 
approach for the effective dynamics. The figures are trivially 
rescaled under the scaling. 



small by the reference of the matter momentum. It is 
then attempting to regard not only <f> as the "internal 
clock" (emergent time) but also p^ as the "internal rod" 
— namely, the measurement of both temporal and spa- 
tial geometries makes sense only in the presence of matter 
content. This observation may support the ideas of the 
relational interpretation of quantum mechanics with real 
rods and clocks such as studied in [l3[ (see also pHII3|). 
although the link is far from clear. If this concept is taken 
seriously, in return, we might be able to further improve 
the /i-scheme to better reflect the underlying physics of 
LQG such that the difference equation of evolution in 
the full LQC theory also respects the scaling invariance 
mentioned above. 
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APPENDIX A: EFFECTIVE DYNAMICS IN 
^-SCHEME 

One of the virtues of the improved strategy (/2-scheme) 
in the isotropic model is to fix the serious drawback in the 
old precursor strategy (^ -scheme) that the critical value 
of the matter density at which the bounce occurs can be 
made arbitrarily small by increasing the momentum p^. 

As the directional densities qi play the same role in 
the Bianchi I model as the matter density p does in the 
isotropic case, we expect that the critical value of gj at 
which the bounce occurs can be made arbitrarily small 
by increasing the momentum p^ in /zo-scheme but is inde- 
pendent of in /2-scheme. The latter is what is shown in 
the main text of this paper. For comparison, the effective 
dynamics in /^ D -scheme is presented here. 

In the effective theory of /x -scheme, we take the pre- 
scription to replace ci by sin(/i°cj) / p,° with the fixed 
numbers fij for discreteness. Analogous to (|3.4p , we have 
the effective (rescaled) Hamiltonian constraint: 



Mo 2 
1 



(Al) 



87rG 7 2 



sin(/^c 2 ) sin(^|c 3 ) 



P2P3 + cyclic terms 



Again, the equations of motion arc given by the Hamil- 
ton's equations and the constraint that the Hamiltonian 
must vanish: 



dp<p 
~dF 

dtf 

IF 



{p^H llo }=0 
W>^„} = P<t>, 



p$ is constant (A2) 
(A3) 



{c u H flo } = 8nG 7 dH ^ 



dpi 



-7 



14 



( S111/JC2 sin (/jgcg) 

X P2 + ~ P3 

V Mi M3 



(A4) 



dpi 
dt> 



{Pi,H„ } 
7~V cos(/z°ci) 



M2 



8.G 7 ^ 



dc± 
sin(/igc 3 ) 



(A5) 



and 



^o( c /.-P/) = ° => P% 
I f sin^ca) sin(^3C 3 ) 



(A6) 



M2M3 



P2P3 + cyclic terms 



From (TAH) and (|A5l). we have 



d_ 

df 



which gives 



sin(/ifc/) 



/'■/ 



= 0, 



sin(/i°cj) 
p/ = 8vrG7ft/C/. 



(A7) 



(A8) 



Taking (|A8|) into (|A6|) gives the same constraints on the 
constant parameters as in (|2. 13[) or (|2. 14[) . 

Substituting ([AS]) into (fA3|) yields 

— % = 8ttGS cos(m?ci)(/C 2 + K 3 ). (A9) 
Pi at' 

By (|A3[) and cos x — — sin 2 x, (|A9|) leads to 



A%=±Vc^g(I^ 
Pi d<p \ k 



2/3' 



1 - 



Qi 



Mo 
0J. crit 



1/2 



• -MO) 

which gives the bouncing solutions with the behaviors 
similar to those given by (|3 . 1 3[) except that the critical 
value of qi at which the big bounce takes place is given 
by 



Mo 
Ql. crit 



-1 3/2 



Kl J p u j 



-, (All) 



which can be made arbitrarily small by increasing the 
value of p$. This tells that /i G -scheme gives wrong semi- 
classical behavior and should be improved by /2-scheme. 
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